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Abstract 



It is established that the distribution of the zero energy eigenfunctions of 
(2 + l)-dimensionaI Dirac electrons in a random gauge potential is described 
by the Liouville model. This model has a line of critical points parameterized 
by the strength of disorder and the scaling dimensions of the inverse partic- 
ipation ratios coincide with the dimensions obtained in the conventional lo- 
calization theory. From this fact we conclude that the renormalization group 
trajectory of the latter theory lies in the vicinity of the line of critical points 
of the Liouville model. 
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Typeset using REVTgX 



Among the problems of localization theory there is one which, until recently, had at- 
tracted less attention than it deserves. This is the problem of spatial correlations of wave 
functions at distances much smaller than the localization length Lc (see Refs. [0,0)- This 
problem is well understood only for extended states, i.e., in the limit of small wave function 
amplitudes t = \ip{x.)\'^. Since extended states explore the entire sample one can neglect 
their spatial variations and treat the Hamiltonian as a random matrix. The distribution 
function P{t) derived by the methods of random matrix theory depends only on the global 
symmetry of the random ensemble and has an approximately Gaussian form (the Porter- 
Thomas distribution; see for example 0). This approach fails for larger fs since the tails 
of the distribution function are determined by rare spatially inhomogeneous configurations 
with high local amplitudes. The first calculation for the asymptotic form of P{t) in small 
two-dimensional samples {L « L^) was performed by Altshuler, Kravtsov and Lerner p| 
using the renormalization group and replicas. The nonperturbative approach based on the 
supersymmetric cr-model had not been used in this context until Muzykantskii and Khmel- 
nitskii |^ pointed out that in order to describe exceptional events most affected by the 
disorder one should look for a saddle point of the supersymmetric cr-model. This idea has 
then been exploited by Fal'ko and Efetov |^ who have derived a reduced cr-model adapted 
to the studies of properties of a single quantum state in the discrete spectrum of a confined 
system and found that the reduced cr-model has a non-trivial vacuum. They have calculated 
P{t) in two loops (the saddle point approximation with Gaussian fluctuations around it). 

Below we discuss only two-dimensional systems. In two-dimensions the conventional lo- 
calization theory of unitary ensemble represents a special case since the first loop localization 
correction to the conductivity vanishes and the localization length is therefore very large: 
Lc ~ exp (Gq), where Gq is the bare conductance (we assume that Go >> !)• Meanwhile the 
behaviour of wave functions becomes nontrivial at much smaller length scales L > exp (Go)- 
In this case the results of Altshuler et al. and Fal'ko and Efetov suggest that the renormal- 
ization flow for P(t), yet eventually turning to the strong coupling (localization), spends 
a lot of "time" in the vicinity of some critical line where the asymptotics of P{t) is given 
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by the log-normal distribution. It is difficult, however, to identify this critical line within 
the replica approach. The supersymmetric saddle point calculations provide a better insight 
since they give the Liouville equation as the cr-model's saddle point condition. 

In this Letter we describe a model of disorder where the line of critical points discovered 
in Refs. |^,^ is stable. This is the theory of (2 + l)-dimensional Dirac fermions in a 
random gauge potential (FRGP model). We show that the distribution function of the 
prelocalized states in this model is described by the Liouville field theory (LFT). This theory 
was introduced by Polyakov ^ in the context of string theory and has been extensively 
studied. We are going to use this accumulated knowledge for the theory of localization. One 
important property of LFT is that it does have a stable line of critical points which seems 
to describe the prelocalized states in the conventional theory of localization. As we shall 
show later, this line is parameterized by the strength of disorder. 

Until present time only two critical disordered systems have been studied: the model 
of a half-filled Landau level (see, for example, ^ and references therein) and the model 
of (2 + l)-dimensional Dirac fermions in a random gauge potential (FRGP model) 



The latter model has an unbounded energy spectrum, but the spectrum of is bounded 
from below {E^ > 0). Since the conductivity of the FRGP model at zero frequency is finite 
and the wave functions with E = have multifractal properties (see Ludwig et al. [H), we 
suggest that the mobility edge of this model coincides with the boundary of the spectrum 
generating operator H^. This conjecture does not generally hold in more complex critical 
disordered systems (e.g. 2D symplectic and ?>D). 

In this letter we study the statistics oi E = Q wave functions of the Abelian FRGP model 
in two dimensions [1^. In the gauge where = t^^dy^^ the Dirac equation has one obvious 



solution with = 0, which we write down in the normalized form: 

— (j'4>(x) 

^.(X) = (1 ± a\^„. ^j^2^.^_2^(x0]l/2 - (1) 

This solution is unique on a closed sphere if the total flux is less than 2. For greater fluxes 
the index theorem predicts the existence of other solutions with E = Q ||T3 |. 



We shall consider V$ as a random variable with the Gaussian distribution: 

P[$] = Zo"^ exp {-^ / d^x[V$(x)]2} . (2) 

Without loss of generality we can choose a = 1 in Eq.(0) and study the moments of the 
distribution function of the first component of the wave function defined by 

iV) = / D$P[$]^?(xi)...^?(x^). (3) 

Notice that all these quantities are invariant under the shift $ ^ $ + c by an arbitrary 
real constant c. The A^-point moment of normalizahle wave function squares can always be 
rewritten as follows: 

G(l.....iV)^ f'^^P^-_°y" /fl«ne--<«.>e-., (4) 

where the action is given by 

S, = jd^x[^^mf + ^^e-^% (5) 

and the exponent exp[— a/i] is introduced to make the integral over /i convergent. Thus, the 
multi-point moment (0) is now expressed in terms of the reducible multi-point correlation 
function of LFT. 

Now we shall recall several facts about LFT that we are going to use. They can be 



found, for example, in Refs. |Tj,|T5|]. To conform to the notations accepted among the field 
theorists, we rescale the field $ = —hcf). We also impose the following boundary conditions 
on the Liouville field: 

0(x) = — Q In |x|^, |x| ^ oo, Q = h+-. (6) 

Thus, we are only considering ground state wave functions that decay algebraically at in- 
finity. These are the prelocalized states. If boundary conditions are not chosen properly, 
the correlation functions of LFT vanish. In the semiclassical saddle point calculations ( 
|l5| JI^JT5[| ) , this feature emerges as a condition for the existence of the saddle point. The 



boundary condition can be formulated as a condition that the total flux through the 
system is equal to Qb. Uniqueness of the wave functions (1) is not violated provided Qb < 2 

There are different ways to implement the boundary condition (^). One can consider the 
LFT on a large disk T of radius R ^ oo and add to the local LFT action (^) a boundary 
term: 

= d'xiWr + 4vr/.e^^1 + ^ /^^ 0d/. (7) 

It is more convenient however to contract the boundary into the point, thereby transforming 
the disk into a sphere. The boundary condition is then equivalent to the insertion of the 
Liouville exponential e"^'^'^'-^^ {R oo) into all correlation functions fl^. For /i = one 



can express correlation functions of the fields e^"*^ in the LFT with boundary conditions in 
terms of correlation functions of the conventional Gaussian model: 



U^) = (8) 



where e is a small parameter introduced for regularization and the subscript Q means that the 
function on the left-hand side is the correlation function calculated with action (|^. Provided 
the neutrality condition = Q is satisfied, all correlation functions are proportional to 

the same factor . We shall see that this factor is compensated by the integral over 

fi. Calculating two-point correlation functions in LFT, we conclude that the conformal 
dimension of the Liouville exponential is 

A(e2°'^) = a(g-a), (9) 



and this is true for /i 7^ also [p!q-|20|| . Note that Q has been chosen such that the exponential 
operator in the Liouville action is marg inal. I.e., A(e2'"^) = 1 thus preserving the criticality 
of the theory. 
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It follows immediately from Eq.(|^) that the conformal dimensions A(g) (g real) of the 
composite operators : ?/'i'^(x) :~ exp[2g60(x)] (the sign : ... : denotes normal ordering) are 
equal to 

A(g) = g (l + ^2 _ 52 \ ^ ^^Q^ 



These dimensions depend on the continuous parameter b which represents the disorder 
strength. The LFT model remains critical for any values of b with the central charge given 
by 

C = l + 6g2 = 1 + 6(6 + 1/6)2. (11) 

However, at 6"^ = M, where M = 1,2, the theory possesses an additional hidden SL(2, 
R) symmetry (see, for example, 0|). It is also known ||15| that three point correlation 
functions have resonances at these values of 6. 

Eq. ( |T0| ) reproduces the dimensions obtained for FRGP in Ref . by the replica trick and 
also the dimensions obtained for the conventional localization theory (in the notation 
of Ref. 1^ 6~2 = Svr^/^z/D). Needless to say that in the conventional localization theory the 
parameter 6 undergoes renormalization towards strong coupling and therefore the described 
equivalence holds only in the crossover regime. As we have mentioned earlier, an extended 
crossover region exists only for the unitary ensemble {jS = 2). 

An important application of (p!OD is to the calculation of the scaling with respect to R of 
the inverse participation ratios: 

(jj^d\^i^(x)) oc R'^^'^\ T{q) = 2(1 - b\){q - 1). (12) 



The applicability of this scaling for large q has been discussed in Refs. ^ and |]T0| , |2T| . It 
was established that r(g) must be a monotonously increasing function of q. This means that 
Eq.(p!2D is at most valid for q < {1 + b^'^)/2. This condition selects operators exp(2g60) with 
qb < (6+ l/6)/2 = Q/2. Let us note that in the weak disorder limit 6 << 1, these exponents 
can still be used for the description of very high participation ratios. 
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This selection resembles the following well known fact from the LFT. Namely, there is 
a difference between operators exp(2a0) with a < Q/2 and a > Q/2 |T^: only the first 
ones (so-called microscopic operators) correspond to local states whereas the latter (so- 
called macroscopic operators) create finite holes on the random surface. In LFT the field 
exp(260) is interpreted as a metric of a two-dimensional surface. In the semiclassical limit 
h « 1 correlation functions of microscopic operators can be obtained by the saddle point 
approximation. The solution of classical equations of motion describes a surface spanned 
on an disk with radius R oo. This surface has constant negative curvature metric with 
spikes (integrable power law singularities of the metric field Qzz ~ kl"''; r/ < 2) at the 
points of insertion of operators with a < Q/2 (for a = Q/2 we have punctures, i.e., 
singularities) . For macroscopic operators there is no classical solution corresponding to a 
surface with a single boundary at infinity which means that each insertion creates a hole of 
a finite size. Since the exponents with q > {1 + b~'^)/2 correspond to macroscopic operators, 
it is not correct, in our opinion, to use them for description of higher order participation 
ratios. It is important to mention, however, that macroscopic operators appear in fusion of 
microscopic ones. 

The authors of Ref. [^] have stressed the importance of the non-universal prefactor on the 
right-hand side of (|l^) to assess the self-consistency of the scaling analysis. This prefactor 
is controlled by the dependency on /i of the correlation functions of the LFT. Let us note 
that before we discussed the results which were independent of /i, so one can put /x = 0. 
Now we have to study the effects which arise in the full theory with nonzero /i. 

It is easy to find the scale (/x) dependence of any correlation function in Liouville theory 

mm- 

G(l,...,iV) = (ne2"'<^(-»))Q = 

i=l 

(vr/i)W-Sti°-)/''F,^...„^(xi, ■ ■ ■ ,xjv). (13) 

Here Fq,j...q^(xi, xat) must be calculated in the Liouville theory with /i = 1 and the 
answer for the generic correlation function is unknown. For three- and four-point functions 
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the answer was found in |2^JT^ for any b. However for the special case when the factor {Q — 
is an integer, i.e., when the correlation function (|T3|) is proportional to an integer 
power of fi, all multi-point correlation functions can be obtained explicitly by expanding 
exp(-S'^,,Q) [cf. (0)] in powers of /i. 

Thus, unlike the case of conformal dimensions where Eq.(^) holds for any b, the multi- 
point correlation functions can be easily obtained only for discrete values of the disorder 
strength 6^ = 1/M, M = 1,2,..., for which the expansion over /i contains only one non- 
vanishing term. Indeed when b~^ is an integer, the multi-point correlation functions of 
the Liouville theory can be expressed in terms of correlation functions of the free bosonic 
field (the Dotsenko-Fateev construction ||23|). To see this, we recall that for the correlation 
function (|^) all ctj = b. Hence, 

N 

{Q-J2ai)/b=l~N + l/b''. (14) 

i=l 

Provided = M with integer M > N — 1 the neutrality condition is fulfilled with the term 
~ ^^^^^^ . The result for the A^-point function on the infinite plane follows immediately: 

,.l-N+M N 1 „ n-N+M \ 

/l-N+M N -I \ /l-N+M 1 \ 

X n n^^xpM n as) 

\ 1=1 i=l K« si I / \ m<n ISmnI / 

(here z.,z = x±i?/ are complex coordinates). The expression for a finite system can be 
obtained by conformal transformation according to the general rules of conformal theory. 
Substituting Eq.(|l5D into Eq.(^) we find that all correlation functions are proportional to the 
same nonuniversal factor D = e^^R'"^^^ d/i/i^/''^e~°'^. We remark that by regularizing 
the theory with the introduction of exp[— a/i], we overestimate the normalization factors 
of the wave functions. In other words, we restrict the allowed spatial fluctuations of the 
wave function amplitude from above, and thereby eliminate the contributions from localized 



states to the statistical average in (^. Finally we note that as shown in Ref. ||2T|, 6 = 1 is 
the largest value of the disorder strength for which quenched and annealed averages [with 
respect to the Gaussian distribution (2)] of the normalization / d^xexp[— 2$(x)] aggree. On 
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the other hand, 6 = 1 corresponds to the minimum of the Liouville central charge Cmin = 25. 
We recall that when 6=1 the total flux through the system reaches the value of 2 which 
corresponds to the (first) change in the ground state degeneracy. 

We conclude with a remark about a difference between even and odd M in Eg. ([T5|) . For 
even M, the operator with maximal dimension A^ax = 1/2 + M/4 is 0^/2(2;) = exp(M60). 
/^From Eq. ([l5D we find: 

ln(|^i2|/a). (16) 



Thus, the two-point correlation function of this operator contains a logarithm. This does not 
happen when M is odd. Indeed, the operator with maximal dimension Amax = (M+1)^/4M 
is 0(Af+i)/2(a;) and its two-point correlation function has the conventional form. It may be 
that, as in other theories with logarithms (see Gurarie |^ and Caux et al. |]lll), Eq. (piBf) 



indicates the presence of an additional operator C{x) in the operator algebra of the theory 
such that 

4D 

{C{z,)Om/2{z2)) = - , iM+2 ' 

\Zl2\ 

{C{z,)C{z2)) = 0. (17) 

We are not certain however that the same interpretation is valid for the Liouville theory 
where there is no one-to-one correspondence between operators and states and one therefore 
must be careful formulating the operator expansion. It is also interesting to note that for 
b = M = 1 we have the Liouville theory with central charge Cmin = 25 which is known 
to have logarithmic operators with scaling dimensions one [^]. As we discussed earlier, 
when 6=1 the total flux through the system reaches the value of 2 which corresponds 
to the appearence of a second ground state in our model. The relationship between the 
existence of marginal logarithmic operator and the change of the ground state degeneracy 
is an interesting open question. 



9 



In more complicated models where the disorder depends on several fields we expect 
that the corresponding universality classes will be connected with VF„-gravities. Another 
interesting problem is to find what deformation of the Liouville theory leads to localization 
and thus to reproduce the renormalization group equations obtained by Altshuler et al. [Q. 

We are grateful to J. Chalker, A. Georges, M. Gunn, K. Efetov, M. P. A. Fisher, D. 
Khveshchenko, V. Kravtsov, D. Khmelnitskii, I. Lerner, A. W. W. Ludwig, A. A. Nersesyan 
and especially to V. Fal'ko and E. Fradkin for illuminating discussions, constructive criticism 
and interest for the work. CM. acknowledges financial support from the Swiss National 
funds. 
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